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ABSTRACT 

We present spectroscopy of 15 star- forming BzK galaxies (sBzKs) with Kab 23 in the 
Subaru Deep Field, for which Ha and some other emission lines are detected in 0.9 to 2.3 /im 
spectra with a resolution of i?=500. Using Ha luminosities, we obtain star formation rates 
(SFRs), and then specific SFRs (SSFRs) dividing SFRs by stellar masses, which are derived from 
SED fitting to BVRi'z'K photometry. It is found that sBzKs with higher stellar masses have 
larger SFRs. A negative correlation is seen between stellar mass and SSFR, which is consistent 
with the previous results for z ~ 2 galaxies. This implies that a larger growth of stellar mass 
occurs in less massive galaxies. In addition, gas-phase oxygen abundances, 12+log(0/H), are 
derived from the ratio of [N II](A6584) to Ha using the N2 index method. We have found a 
correlation between stellar mass and oxygen abundance in the sense that more massive sBzKs 
tend to be more metal rich, which is qualitatively consistent with the relation for UV-selected 
z ~ 2 galaxies. However, the metallicity of the sBzKs is ~ 0.2 dex higher than that of UV-selected 
galaxies with similar stellar masses, which is significant considering the small uncertainties. The 
sBzKs in our sample have redder R — K colors than the UV-selected galaxies. This galaxy color- 
dependence in the oxygen abundance may be caused by older or dustier galaxies having higher 
metallicities at z ~ 2. 
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1. Introduction 

Recent studies suggest that the era of z ~ 2 is a 
turning point in galaxy formation and evolution. 
The cosmic star formation rate (SFR) begins to 
drop at z ~ 1 — 2 from a flat pl a teau at higher 



redshifts ( Dickinson et al. 20031 : iFontana et al 



20031 iLv et all 120071 ). Significant evolution of 
the Hubble sequence occurred at z ~ 1 — 2 



(jKaiisawa fc Yamadal l2001'). It is also found that 
the number density o f QSOs has a peak at z ~ 2 



( Richards et al. 20061 ). These facts suggest that 



dramatic changes in the galaxy population oc- 
curred at z ~ 2, which are important to study in 
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further detail. 

Galaxy properties are determined from spec- 
tral characteristics such as emission lines, ab- 
sorption lines and strong continuum breaks. At 
z ~ 2, the strongest of these features fall in the 
near-ultraviolet (NUV) and near-infrared (NIR) 
wavebands. NUV and NIR observations are thus 
essential to reveal detailed properties of z ~ 2 
galaxies, but these observations encounter many 
difficulties, including poor sky transparency ex- 
cept for some atmospheric windows, and strong 
OH airglow emission lines and thermal emission 
from the atmosphere. This is why the era of 
z ~ 2 had been called a 'redshift desert' until 
recently. However, recent advances in wide-field 
imaging in the NUV and NIR, along with multi- 
object spectrographs, are beginning to fill in this 
'desert' (e.g., Subaru Deep Field (SDF) survey, 
GOODS-Nort h and South, UKIDSS, GMASS . 
and MUSYC dMcCarthv et alJll999t iReddv et al 



Z correlation at z ~ 2 as there is in the local uni- 
verse, with the more massive galaxies being more 
metal rich. Although the correlation has a similar 
slope to the local one, it is shifted to lower metal- 
licity by ~ 0.3 dex. They have also found that 
SFRs derived from Ha luminosities are in good 
agreement with those from de-reddened rest-UV 
luminosities, with an average dust-corrected SFR 
of - 3OM yr- x . 

However, the rest-UV color selection misses 
about half of z ~ 2 galax ies, which have a large 
amount of dust extinction ( Kong et al . 2006). The 
spectroscopic properties of this dusty population 
have yet not been clearly revealed, if-limited sam- 
ples are better suited than UV-limited ones to in- 
vestigate the relations between stellar mass and 
other fundamental properties at z ~ 2, because 
they approximately correspond to samples limited 
by stellar mass, witho ut excluding galax ies with 



2006b; Erb et al. 2006c; Hayashi et al. 200 7tlDaddi et all 
2007tlGrazian et al.l2007tlLane et alj2007t lHalli dav et al 



200RtlKriek et al.ll2008L ILv et alJl2008h ). 

A photometric selection technique using B — z 
and z — K colors has been proposed to choose an 
unbi ased sample of z ~ 2 galaxies (BzK galax- 



ies 



Daddi et al 



(2004)). Using these colors to 
find galaxies with a Balmer break between z and 
K, the selected galaxies are then classified into 
star forming galaxies (sBzKs) and passively evolv- 
ing galaxies (pBzKs). Thanks to wide and deep 
imaging surveys in the NIR as well as the optical, 
large new BzK samples are enabling us to mea- 
sure their statistical propert i es, such as clustering 
strength dKong et all l2006t lHavashi et all 12007 : 
Blanc et al.l 2008). The multiwavelength data 
also show that sBzKs are very actively formin g 
stars! Dannerbauer et al. 20061 : Daddi et al. 2007 ). 
On the other hand, due to lack of spectroscopic 
observations, there are relatively few accurately 
known redshifts, which are needed to derive accu- 
rate stellar masses and star formation rates from 
SED fitting to the multiwavelength data. Also, 
their physical properties such as metallicity are 
still not well-known. 

More spectroscopy has recently become avail- 
able for z ~ 2 gala xies selected by their UV colors. 
lErb" et al. (|2006al Tbh presented a mass-metallicity 
(M-Z) relation and SFRs for UV-selected galaxies 
at z ~ 2 using 114 NIR spectra. They found a M- 



large dust extinction. iKriek et alj (|2008f ) carried 



out a NIR spectroscopic survey for if-bright galax- 
ies at z ~ 2.3. They suggest that studies with 
broadband photometric data alone may overesti- 
mate the number of massive galaxies at z ~ 2, and 
underestimate the evolution of the stellar mass 
density, again indicating the importance of spec- 
troscopy. 

In this paper, we present NIR spectroscopy for 
44 BzK galaxies in the SDF, and their spectro- 
scopic properties. The structure of this paper is as 
follows. The photometric and spectroscopic data 
are described in Sj2j In $3l we examine the proper- 
ties of the detected emission lines, including spec- 
troscopic redshifts and fluxes. We then investigate 
active galactic nuclei (AGN) contamination in our 
spectroscopic sample in 21 In SJSJ physical proper- 
ties are derived from the emission lines, and then 
discussed. A summary is given in Sj6] Through- 
out this paper, magnitudes are in the AB system, 
and we adopt cosmological parameters of h = 0.7, 
f! m o = 0.3 and f^Ao = 0.7. 

2. Data 

2.1. BzK galaxies in the SDF 

Here we briefly describe the imaging data and 
sample selection of BzKs in the SDF, which 



are al most the same as those in lHavashi et al 
(|2007l ). The SDF is a large "blank" field cen- 
tered on (13 h 24 m 38!9, +27°29'25'/9; J2000), a 
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rectangle of 29.'7 x 36.'7, with deep multi-wave- 
length data from the NUV to mid-IR (MIR) wave- 
bands. Among these data, the optical and the 
NIR data are mainly used in this study. The op- 
tical data, B, V, R, i', z', NB816, and NB921, 
were obtained with the Suprime-Cam on the Sub- 
aru telescope for the S ubaru Deep Field Project 



( Kashikawa et al. 20041 ). and then the NIR data, 



J and K, were obtained with the wide field cam- 
era (WFCAM) on the United Kingdom Infra-Red 
Telescope (UKIRT) (Motohara et al. in prepa- 
ration). The NIR data available to this study 
were only what was obtained on 2005 April 14 
and 15, which limited the overlapping area with 
the optical data to about 40% of the whole SDF 
([Havashi et al. Il2007t> . The 3er limiting magnitudes 
in a 2"-diameter aperture are 28.45 magnitude for 
B, 26.62 for z', and 24.05 for K. The optical im- 
ages have been convolved to a seeing size of 1'.'14, 
which is the value for K. 

MIR data, 3.6 - 8.0 /im, were obtained with 
the InfraRed Array Camera (IRAC) on the Spitzer 
space telescope. Since the point spread functions 
in the MIR data are much broader than those 
in optical and NIR data, some galaxy images are 
blended with nearby objects. For these, we had to 
remove the contaminating flux contributions from 
close objects in order to do accurate photome try. 
To do this we use GALFIT (jPeng et alJl2002h to 
deblend images with adjacent objects (Hayashi et 
al. in preparation), and accurately measure the 
MIR fluxes of objects. The 3er limiting magnitudes 
are 23.60 for 3.6/im, 23.32 for 4.5/im, 21.65 for 
5.8/im, and 20.99 for 8.0/im, which are measured 
with 4.4", 4.4", 5.0", and 6.4" diameter apertures, 
respectively. 

We applied the Bz K color selection proposed 
by iDaddi et all (|2004f ) to the B, z' , and K data 
in the 180 arcmin 2 of the SDF which have good 
detection uniformity in all three passbands. In 
selecting BzK galaxies, deep optical data are cru- 
cial, because objects without B or z' detections 
can take arbitrary positions on the BzK color di- 
agram, with almost all of them left unclassified. 
We selected the BzKs with great c are for the op 



tically undetected objects (See also lHavashi et al 
(|2007l) ). We obtained samples of 1092 sBzKs and 



2.2. Spectroscopic observation 

NIR spectroscopy of the BzKs in the SDF was 
carried out with the Multi-Object InfraRed Cam- 
era a nd Spectrograph (MOIRCS; llchikawa et al 



2006) on the Subaru telescope on 2007 May 3 and 
4. MOIRCS has a wide field of view of 4' x 7', and 
its z J500 and HK 500 grisms provide spectra with 
a resolution of i?=500, covering wavelengths of 0.9 
- 1.8 /im and 1.3 - 2.5 fim, respectively. The dis- 
persions of the grisms are 5.57 and 7.72 A/pixel, 
respectively. 

Forty-four spectroscopic targets were selected 
from the samples of 1092 sBzKs and 56 pBzKs. 
We gave priority to K < 22.0 BzKs in our samples 
for slitmask design, and then K > 22.0 BzKs were 
allocated to the blank remaining slits. We used 
two slitmasks to cover the 44 BzKs, that is, each 
slitmask covers 22 BzKs. Each slit has a width of 
0"7 with a length of 10 - 12 " according to the 
space between the slits. Among the 44 galaxies, 
40 are classified as sBzKs and 4 as pBzKs. Half of 
them are fainter than K = 22.0. Figure Q] shows 
the colors of the 44 BzKs. 

We observed them consecutively at two posi- 
tions (A and B) on the slit, using the zJ500 and 
#1^500 grisms for one slitmask per night with sin- 
gle exposure times of 600 or 900 sec. The total on- 
source integration time was then 2 - 2.6 hours for 
both grisms. During the two nights, the weather 
was good, and the seeing was as good as ~ 074 
in K, except for the end of May 4, when it was 
> 1"0. We also observed AO stars with V ~ 9.0 
magnitude at different air masses with both grisms 
on each night to correct for the telluric absorption 
and the instrumental efficiency. Before and after 
the observations, dome flats with lamps on and 
off, and Th-Ar lamp data were obtained for wave- 
length calibration. 

2.3. Data Reduction 

Data reduction is done with standard iraeEI 
procedures. First, we create A - B frames from 
successive frames observed at the two positions. 
The A - B frames are then divided by a flat- 
field image, and bad pixels and cosmic rays are re- 



56 pBzKs with K < 23.2, corresponding to the 5a 
limiting total magnitude. 



9 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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Fig. 1. — The colors of 44 targets. The filled cir- 
cles are BzKs with K < 22.0, and the filled trian- 
gles are those with K > 22.0. BzKs with emission 
lines detected are squares. 



moved. The flat image is created by taking differ- 
ence between the dome flat frames with the lamps 
on and off, and normalizing it. Next, distortion is 
corrected using calibration data files provided by 
the MOIRCS instrument team. Wavelength cal- 
ibration is done with only the OH airglow lines 
for the HK data, since there are numbers of the 
strong OH lines in the wavelengths longer than 
1.1/xm. On the other hand, for the zJ data, Th- 
Ar data are also used for wavelength calibration at 
wavelengths shorter than 1.1/xm. The wavelength 
uncertainties in the strong lines are less than 15A. 

Residual sky subtraction is then carried out, 
since only the A - B procedure may not com- 
pletely remove the sky background due to its time 
variation. The telluric absorption and the instru- 
mental efficiency are corrected using the spectra 
of AO star s and the stellar spectral library given 
bv lPickled (1998). Finally, all the frames for each 
galaxy are added and ID spectra are extracted 
combining 5 pixels along a slit. 

Flux calibration of the spectra is done using 
the photometry from the imaging data in z' and 



spectrum. 

3. Measurement of emission lines 
3.1. Detection of emission lines 

We have detected Ha and some other emis- 
sion lines in 15 of the 40 sBzKs, 9 of which have 
K > 22.0. Since there are many strong OH air- 
glow lines in the NIR region, we have first con- 
firmed all lines in the 2D spectra by visual inspec- 
tion. Among the confirmed lines, then, we regard 
lines with a peak flux larger than the 3<r sky noise 
as a detection using ID spectra with a resolution 
of i?=500, that is, unbinned spectra. For the 4 
pBzKs, no spectra with enough S/N to detect any 
Balmer/4000A break were obtained due to the in- 
sufficient exposure times. In what follows, we ex- 
amine and discuss spectroscopic properties of the 
15 sBzKs with emission lines detected. Figures [2] 
- 0] show the binned ID spectra of the lines, and 
Tables [1] and [5] show properties of the 15 sBzKs. 

To measure spectroscopic redshifts and fluxes 
of the emission lines, Gaussian profiles with a con- 
stant continuum are fitted, where the free param- 
eters are amplitude, line width, redshift, and con- 
tinuum level. The ler error of the spectrum is 
used as a sigma in the calculation of x 2 . Because 
Ha (A6563) and [N II] (AA6548,6584) are blended, 
they are simultaneously fitted, assuming the same 
line width and redshift, and a 6584/6548 ratio of 3- 
to-1. The la error of each parameter in the fitting 
is estimated as the range where the increase in x 2 
from the best-fit value is less than the appropriate 
value for the number of deg rees of freedom (s ee 
§14.5 of Numerical recipes bv lPress et al. ( 1995h ) . 



K. The lcr error of the spectrum is estimated from 
the standard deviation in the sky region of the 2D 



3.2. Spectroscopic redshifts 

Spectroscopic redshifts are derived from the Ha 
lines using the central wavelength of the fitted 
Gaussian. For four sBzKs, which are ID21360, 
ID21667, ID22602, and ID27677, only one emis- 
sion line is detected, but we regard it as Ha for the 
following reasons: (1) The redshift of ID21360 is 
also confirmed to be 1.174 by optical spectroscopy 
with t he Hectospec o n the Multiple Mirror Tele- 
scope (|Lv et al.l 120081) . Therefore, the line is cer- 
tainly Ha. (2) For the other three objects, it 
is possible that the detected lines, seen in the 
HK grism spectra, are from transitions other than 
Ha. If the lines are Ha, the redshifts should be 
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Fig. 2. — The spectra of Ha and [N II] emission lines we have detected. These spectra have a binned 
dispersion of ~ 16 A. The solid lines show the spectra and the shaded regions in gray are the la sky noise. 
We regard lines with a peak flux larger than the 3<r sky noise as a detection using unbinned ID spectra. 
Labeled redshifts are derived from Ha lines. 



1.75, 1.87, and 1.71, respectively. On the other 
hand, if the lines are [O III](A5007), those would 
be 2.60, 2.76, and 2.55, respectively. The transi- 
tion of H/9 would result in even higher redshifts. 
Then, because the BzK technique isolates select- 
ing star forming galaxies at 1.4 < z < 2.5 by 
detecting the Balmer break between the z-band 
and K, Ha is the most likely identification of lines 
in the HK region. (3) In addition, we note that 
almost all the derived spectroscopic redshifts are 
consistent with photometric redshifts. The photo- 
metric redshifts are obtained by fitting template 
SEDs to BVRi'z'K phot ometric data using th e 
Hyperz code (version 1.1) (jBolzonella et al.l l2000) . 
The used set of template SEDs is what is pro- 
vided in the Hyperz package, which consists of the 
observed spectra bas ed on the CWW templates 
(jColeman et al.lll980h . The standard deviation of 



&{z sp ec — -Zphot) is 0.61 for the 12 sBzKs with the 
secure spectroscopic redshfits. The photometric 
redshifts for ID21667 and ID22602 are 1.47 and 
1.46, respectively, supporting our identification of 
their lines as Ha. ID27677 would have a photo- 
metric redshift of 2.46, suggesting the possibility 
that the line is H/3 or [O III]. Since the flux ra- 
tio of JO HI] to H/3 can widely r ange from 0.1 
to 10 ( Veilleux fc Osterbrockl 1987 ). it is possible 
that only one of the two lines is detected. Note 
that [O III](A4959) would not be detected due 
to the inadequate S/N even if the identification 
of [O III](A5007) is correct. These considerations 
suggest that all the three of Ha, H/3, and [O III] 
are candidates for the single line of ID27677, and 
that no secure conclusion is available for the line 
identification. Therefore, adopting the considera- 
tion (2), we assume that the line for ID27677 is 
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also most likely to be Ha. Even if the line is not 
Ha, our conclusions discussed below would remain 
unchanged. 

Since the spectral resolution is i?=500, uncer- 
tainties in spectroscopic redshifts are estimated to 
be Az = AA obs (Ha)/6563 - 0.002. 

Figure [5] shows the redshift distribution of the 
sBzKs with spectroscopic redshift, indicating a 
somewhat unexpected feature. The distribution 
does not peak at z ~ 2, and no object lies at 
z > 2, contrary to the fact that BzK color selec- 
tion is proposed to select g alaxies at 1.4 < z < 2.5. 
In fact. lDaddi et al.l (|2007l ) have reported that the 
redshift distributions of sBzKs with spectroscopic 
and photometric redshifts both peak at z ~ 2, and 
that this selection technique works well to select 
z ~ 2 galaxies, which is also supp orted by other 



resul t s with photometri c redshifts (jGrazian et al 



4700 4800 4900 5000 5100 



20071: iBlanc et all 120081 ) . The discrepancy of our 
result may be due to a flux- limited bias, or the 
fact that the strongest emission line, Ha, is shifted 
into a wavelength region of high sky background 
for z > 2, as discussed in the next section. The 
opposite problem also occurs, namely that some 
galaxies at z < 1.4 are selected, which are also 
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Fig. 5. — The redshift distribution of the sBzKs 
with spectroscopic redshift derived from the Ha 
lines. Errors are Az ~ 0.002. 



seen in Fig. 2 of iDaddi et all (|2007l) . These lower- 
z objects are near the boundary of BzK selection 
criterion, and lie 2.5 - 2.7 sigma from the color 
selection threshold. This result suggests that a 
fraction of star-forming galaxies at lower redshifts 
(z < 1.4) can meet the sBzK criterion. 

3.3. Objects without detected emission 
lines 

Among the sBzK targets, no emission line is de- 
tected for 25 sBzKs. It is important to know the 
cause of the non-detections. Judging from the dis- 
tribution of these sBzKs without a line detection 
on Figure [TJ it is possible that they ha ve large dust 
extinc tion. Applying equation (4) of IDaddi et al. 
(2004) to the B — z color, we find that sBzKs with 
E(B — V) larger than ~ 0.5-which corresponds to 
Av=2.0 magnitude-have no lines detected. This 
may be one of the reasons for the non-detection 
of lines. Since there are many strong OH airglow 
lines in the NIR region, it is also possible that an 
emission line falls on a strong OH line. In addition, 
lines are redshifted into the if-band region in the 
case of z > 2 sBzKs. As described in section I3~2l 
the S/N in the region may be too low for lines to 
be detected. In fact, there were no lines detected 
at a wavelength longer than 2.0/xm. Finally, it is 
also possible that only sBzKs with higher SFRs 
have strong enough lines to be detected in this K- 
limited sample. The lowest Ha flux we detected is 
4.2 xlO -17 ergs s _1 cm -2 . 
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4. AGN contamination 

Before considering the spectroscopic properties 
of sBzKs, we should account for possible AGN 
contamination in our sample. We use two meth- 
ods to identify AGN candidates. One is the BPT 



emiss ion-line diagnostic diagram (Ba ldwin et al 



19811 ). and the other is a power-law shaped SED 



from optical to MIR wavelengths. 

Figure [6] is the emission-line diagnostic diagram 
showing the [O III] (A5007)/H/3 vs. [N II] (A6584)/Ha 
ratios. The only three sBzKs-ID26821, ID26840, 
and ID31569-have all the required lines detected. 
Undetected line fluxes for the other galaxies 
are estimated as follows. We derive H/3 fluxes 
from the Ha fluxes, E(B — V) described in §5.1, 
and the assumed intrinsic Balmer decrement of 
Ha/H/3 = 2.86. The fluxes of [O III] and [N II] are 
the ltr upper limits derived from the sky noise at 
the wavelengths where each line should be seen, 
assuming the same line width a s Ha As dis- 
cussed in iKewley fc Dopital (|2002l ) and lErb et al.l 



(|2006al ). the flux ratio of [N II] /Ha is enhanced 
by shock excitation (in LINERs), or by the hard 
continuum which photo-ionizes the Narrow Line 
Region (NLR) of all Se yfert nuclei. Figure 7 in 
Kewlcy fc Dopital (|2002h shows that the ratio of 
[N II] /Ha is always less than 0.63 in purely star- 
forming galaxies. Among the 15 spectroscopically 
confirmed sBzKs in our sample, four- ID21193, 
ID22424, ID25528, and ID31340-have [N II] /Ha 
greater than 0.63, which may imply a significant 
contribution from either a LINER or Seyfert NLR 
component. In particular, ID25528 shows the high 
[N II] /Ha ratio of 1.2, suggesting that it is an 
AGN. On the other hand, we cannot strongly as- 
sert that the other three galaxies are LINERs or 
Seyferts due to the uncertainties of their line ra- 
tios. Also, as described below, the SEDs suggest 
that their continuum is dominated by stellar emis- 
sion, and that the AGN contribution is small. If 
they are type II AGNs, the strongest observed 
contribution to the spectrum from the active nu- 
cleus would be high-ionization emission lines, not 
continuum. Similarly, one galaxy-ID 31569-shows 
a high enough [O III](A5007)/H/3 ratio of 4.9 to 
be a possible Seyfert 2 galaxy. However, the weak 
[N II] line suggests that this is instead a metal- 
poor star-forming galaxy. 

These ratios of forbidden to Balmer emission 



lines do not however provide a definitive search 
for AGN. This is because a strong contribution 
from the Broad Line Region (BLR) in Seyfert 1 
nuclei produces a low [N II] /Ha ratio which would 
not be noticeably different from those of H II 
regions. Our spectral resolution of 600 km s -1 
should be sufficient to detect most broad wing of 
Ha, but no such wing is found. Furthermore, only 
one o f the Ha lines-ID 26821-satisfies the crite- 
ria of iHicks et all (|2002h set for an Seyfert 1: EW 
> 100A and log L(Ha, observed) [ergs s _1 ] > 42.5 
(Table[2|). This implies that there is a contribution 
from a Seyfert BLR component in ID26821 sBzK. 
However, as shown in fl5.2.2[ the SFR derived from 
Ha luminosity is in agreement with that from UV 
luminosity within a factor of 2, which may indi- 
cate that the contribution from a BLR to the Ha 
line is not significant. 

We have also applied another test for Seyfert 
1 nuclei in our sample. We checked the long- 
wavelength SEDs of the 15 sBzKs using the pho- 
tometry from 3.6 — 8.0/j.m taken with IRAC on 
the Spitzer space telescope. Most of the sBzKs 
have a strong bump at a rest wavelength of 1.6 
yitm in their SEDs, indicating th at their SEDs 
are dominated by ste l lar em ission ( Sawicki 2002; 
lAlonso-Herrero et ail l2007t ). AGN-dominated 



galaxies, especially those with broad emission lines 
(Seyfert l's) have a power-l aw SED at rest wave- 
lengt h s longer than 1 um (Malkan fc Filippenkol 



Il983t ISpinoglio et~aH 119951 1 In our sample, 
no such clear power-law SEDs are seen, but 
three sBzKs have power-law-like SEDs: ID21973, 
ID26821, and ID27677 (Figure El). These might 
perhaps contain "Seyfert 1.9" or "Narrow Line 
Seyfert 1" nuclei. However, intense star formation 
also produces substantial near-IR emission from 
heated dust, which could mimic such an SED. Fig- 
ure [7] shows these three power- law-like SEDs with 
SED models of st arbursts overplotted (SBURT: 
Takagi et aDl2003h . These model SEDs have no 



contribution of AGN. 

With all these considerations, we classify 
ID25528 and ID26821 as AGN. Although the AGN 
contribution to the spectra, particularly the con- 
tinuum, may not be significant, the two galax- 
ies are excluded from the discussions on emission 
lines in the following section. We then classify the 
three sBzKs of ID21193, ID22424, and ID31340 
as AGN candidates judging from the [N II] /Ha 
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Fig. 6. — The emission-line diagnostic diagram showing the [0 III] (A50Q7) /H/3 vs. [N II](A6584)/Ha ra- 
tios. The solid line is the boundary separating star forming galaxies (left side) and AGNs (right side) 
(Kauffmann ct al. 2003). Circles show objects identified as star forming galaxies, while a square and trian- 
gles as an AGN and possible AGN candidates, respectively. The undetected fluxes of [O III] and [N II] are 
estimated from the la upper limits, which are shown by arrows. The Hf3 fluxes are inferred from the Ha 
flux and E(B — V) on the assumption of an intrinsic Balmer decrement of Ha/H/3 = 2.86. The error bars 
are shown with la. 



ratios. However, since it is possible that the three 
AGN candidates are star forming galaxies, we 
consider them as well as the other ten sBzKs clas- 
sified as star forming galaxies in the following 
discussions. 

( 20061) 



Kong et al 



have estimated that about 
25% of sBzKs are AGNs using XMM-Newton X- 
ray data, which would imply that 3-4 sBzKs in 
our sample c ould b e AGN. However, the sample of 
Kong et all (120061) has only if -bright BzKs. Ac- 
cording to I Alonso-Herrero et al. I (120071 ). the char- 
acteristic stellar mass of the host galaxy of AGN 
increases with redshift, and at z ~ 2 it is ~ 
(1 — 3) x 10 11 M Q . In our sample, almost all galax- 
ies are less than 10 n MQ, which indicates that the 
fraction of AGN c an be lower than 25% . The NIC- 
MOS imaging of IColbert et al I (|2005l ) also indi- 



cates that fewer than a quarter of these less mas- 
sive galaxies should host an AGN detectable in 
X-rays. The AGN fraction of 2/15 (13%) in our 
sample is consistent with the number of AGN ex- 
pected from other studies. If the three candidates 



are actually AGNs, the fraction increases to 33%, 
which would still be consistent with other studies. 

5. Properties of star-forming BzK galaxies 

Emission line information enables us to esti- 
mate SFR and metallicity. In the local universe, 
more massive galaxies tend to have lower SFR s 
and higher metallicities (e.g. JGallazzi et aLll2005h . 
It is important to examine whether a similar rela- 
tion exists among these properties for high redshift 
galaxies. 

In this section, we first describe the properties 
obtained from multiwavelength photometric data 
by SED fitting method. Then, spectroscopic prop- 
erties are discussed, along with the photometric 
properties. 

5.1. Properties from SED fitting 

Properties such as stellar mass and dust ex- 
tinction, E(B — V), are obtained from fitting 
synthetic model spectra of continuous bursts 
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Fig. 7. — The long-wavelength power-law-likc 
SEDs of ID21973, ID26821, and ID27677 sBzKs, 
where model SEDs of starbursts (w i thout any 
AGN contribution) from iTakagi et aD ( 2003 ) are 
overplotted. 



(|Bruzual fc Chariot |[200^) to the BVR i 'z'K pho 
tometry, assuming a standard ISalpeterl (|l955l ) ini- 
tial mass function (IMF), solar a bundance, and 



the e mpirical dust extinction of iCalzetti et al 
(2000). These are shown in Table [TJ Stellar mass 



is one of the most important parameters of galaxy 
properties and the most robust one obtained from 
SED fitting. Although sBzKs are well-fitted by 
the continuous burst model, fits using the r model, 
where the SFR declines exponentially with time, 
would make the stellar masses slightly smaller 
than those with the continuous burst model. Note 
that SEDs of the simple stellar population (SSP) 
model are irrelevant in this work, since sBzKs are 
still star- form i ng at z ~ 2. It is also noted that 
if a IChabrierl (|2003h IMF is used, stellar masses 
(and SFRs) are smaller by a factor of 1.8 than 
those with a Salpeter IMF. Thus, stellar masses 
are determined within a factor of ~ 2.5. Fitted 
E(B — V)s are mainly used to calibrate the dust 
extinction of emission line fluxes in section 15.21 
As discussed in section 0] almost all SEDs of the 
sBzKs are dominated by stellar emission with a 
clear 1.6/im bump. Therefore, any AGN contribu- 



tion to the SED is small enough that parameters 
obtained from the SED fitting reflect properties of 
the stellar populations. 

5.2. Star Formation Rates 

5.2.1. SFR derived from Ha and UV continuum 
SFRs for sBzKs are derived from the luminos- 



i ty of Ha using the equation given in iKennicutt 
(Il998h : 



SFRHcjMoyr- 1 ] = 7.9 x 10- 42 L H a[ergs s" 1 ], 

(1) 

assuming solar abundance, Salpeter IMF with 
mass limits of 0.1 and 100 M@, and a constant 
SFR. Ha luminosities are corrected for the dus t 
extinction by equation (3) in lCalzetti et a l. (2000), 
E(B-V) = 0ME gas (B-V). The amount of neb- 
ular extinction is then estimated fro m E V!>R (B — V) 



using the extinction law given in ICalzetti et al 
(2000). The median amount of extinction of Ha 
is 2.83 magnitude. 

The amount of extinction at Ha can be directly 
estimated from the apparent Ha/H/3 ratio. How- 
ever, we do not apply this method for our sample, 
because only three galaxies have detection of both 
lines and even for these three the amount of ex- 
tinction is not strongly constrained due to a large 
uncertainty in the H/3 flux. 

The rest-UV luminosity is also used as a tracer 
of SFR. We derive SFR from the continuum at 
1500A, L 1500, which is predicted by the best- fit 
model to the BVRi'z'K photometry. L1500 is cor- 
rected for the dust extinction with Calzetti ex- 
tinction law and converted into a S FR using the 
equation given in lDaddi et al.l ( 2004f k 



SFR uv [Af yr- 



£1500 [ergs s 1 Hz l ] 27 

8^5 Xl ° ■ 

(2) 



5.2.2. Comparison between SFRh q and SFRuv 

Figure [8] compares the SFR from the rest-UV 
luminosity with that from Ha. A correlation is 
seen between the two SFRs, but SFRh q tends 
to be larger than SFRuv by a factor of three. 
The median SFR Hq and SFRuv are 182 and 61 
M yr _1 , respectively. Also, this discrepancy be- 
tween SFR measured by Ha and from shorter 
wavelength observations is quite similar to what 



9 



Hicks et al.l (|2002f) found in their sample of IR- 
selected galaxies. 

On the other hand, Erb et al. ( 2006bl ) report 
that SFRhq for UV-selected z ~ 2 galaxies is in 
good agreement with SFRuv, which is inconsis- 
tent with our result. However, their observed Ha 
luminosities are similar to ours. That is, the me- 
dian Ha l uminos ities of our sample and that of 
Erb et alj (|2006bl ) are 2.0 and 2.2 times 10 42 ergs 
s , respec tively. Th i s disagr eement on resulting 
SFRs with lErb et alj (|2006bl ) may be due to dif- 
ferent corrections for dust extinction at Ha. They 
used the color excess for the stellar continuum to 
de-redden Ha, that is, E gas {B -V) = E(B - V"), 
which implies that their amount of dust correc- 
tion for Ha luminosity is systematically lower than 
ours by 1.6 magnitude on average. This leads to a 
difference in derived SF Rs of a factor of 4 .4. Us- 
ing the same method as lErb et all (l200fibl ) for de- 
reddening Ha, our correlation between SFRh q and 
SFRtjv becomes consistent with that of lErb et al. 
(|2006br ). In this case, the median SFRh q is 50 



1000 



M Q yr- 



Daddi et alJ (|2004h report that most color ex- 



cesses of 24 sBzKs in the K20/GOODS area, which 
are derived from U to K SEDs, fall within a range 
of 0.2 < E(B - V) < 0.5, which is consistent with 
our result (Table [lj. The median E(B — V) in our 
sample is 0.38. This implies that our estimation 
of E(B - V) by SED fitting is reliable. 

This difference between SFRh q and SFRuv im- 
plies two possibilities. One is t hat a Salpeter 
IMF i s unsuitable, if the equation of lCalzetti et al. 
(|2000l ). E(B-V) = 0.UE sas (B-V), is assumed to 



be valid for star- forming galaxies at z ~ 2. A more 
top-heavy IMF than Salpeter would account for 
the difference, since nebular emission results from 
ionizing photons shortward of the Lyman limit 



The o ther is that the equation of ICalzetti et al 
(2000) cannot be applied to z ~ 2 galaxies. In 



any case, in what follows we use SFRh q for sBzKs 
in the comparison with their stellar masses. Even 
if SFRuv were used instead, the discussion in the 
next section would remain unchanged. 

5.2.3. Relation to stellar mass 

We examine the relation between SFRh q and 
stellar mass. Figure [9] suggests that sBzKs with 
higher stellar masses have larger SFRs, which is 




10 100 1000 

SFR corr (rest-UV) [M. yr"'] 

Fig. 8. — Comparison between SFRh q and 
SFRuv- SFRhqS are derived from Ha luminosi- 
ties, and then SFRuvs are from the luminosities 
at 1500 A estimated from the best-fitted SEDs. 
Both SFRs are corrected for dust extinction. The 
dotted line shows that SFRuv is equal to SFRjja- 



qualitatively consistent with iDaddi et al.l (|2007T ). 



However, the slope of our relation is shallower 
than theirs. One cause for this disagreement in 
the slope may be a difference between SFRhq and 
SFRuv ■ Figure [8] suggests the difference between 
these two estimates of SFR is larger at lower SFRs. 
Since less massive sBzKs have lower SFRs, this 
difference results in shallower slope of the correla- 
tion between stellar mass and SFR Hq ■ Other cause 
may be a differ ence in the method used to derive 
stellar masses. IDaddi et al] (|2007| ) have used the 
empirica l relation between t he stellar mass and K 
(|2004l ). which is calibrated 
K color. We compare the stellar 



given m 
based on z 



Daddi et al 



masses derived from a SED-fi t ting w ith those from 



the relation of IDaddi et al.l (|2004t ) , finding that 



the difference between the two estimates of stellar 
mass is larger at smaller stellar mass. The stellar 
masses from the relation are larger by a factor of 2 
- 3 in ~ 10 M . Because this relation is derived 
using bright K < 22 BzKs, the stellar masses of 
faint BzKs may be overestimated compared with 
those from SED fitting, resulting in the shallower 



10 



correlation between stellar mass and SFRh q . 

It is also possible that a flux-limited bias, dis- 
cussed in §§3.2 and 3.3, produces this disagree- 
ment. If galaxies with (E(B - V)} = 0.38 and 
F(Ha) — 4.2 x 10 _17 ergs s _1 cm -2 , which is the 
lowest Ha flux we have detected, are at (z) ~ 1.7, 
the SFR would be ~ 87M Q yr _1 . This suggests 
that it is likely that we cannot detect any emis- 
sion lines for galaxies with SFR < 87 Mqyy" 1 . 
Moreover, it is found that SFRh q s is three times 
larger than SFRuvS on average. Combining this 
with the relation between stellar m ass an d SFRuv 

<l2007h . we find 



1000 



given in Fig. 17 of iDaddi et al 



that the limiting SFR of 87 M Q yr _1 corresponds 
to a stellar mass of ~ 1 x 10 10 M©. These estimates 
indicate that it is possible that no galaxies with 
< 1 x 10 10 M Q and < 87M yr _1 are detected due 
to flux-limited bias. However, even if Figure [9] is 
limited to a region with SFRs higher than the esti- 
mated detection limit, a mild correlation between 
stellar mass and SFR can be seen. If there are 
les s massive ga l axies with lower SFRs, as found 



Daddi et al.l ((2007), the correlation would be 



steeper. 



Erb et al.l (|2006bf ) have also reported a simi- 



lar correlation between SFR and stellar mass for 
UV-selected z ~ 2 galaxies, but their correlation 
is shifted to much lower SFR than ours. This 
is mainly due to differences in the assumed IMF 
as well as dust corre c tion, a s discussed in section 
15X21 In lErb et all (|2006bl ). a Chabrier IMF is 
used to estimate stellar masses and SFRs. SFRs 
derived with a Chabrier IMF are smaller than 
those with a Salpeter IMF by a factor of 1.8. 
These differe nces make o ur SFR 7.9 times larger 



than that of lErb et al.l (|2006bl ). After correct- 
ing these systematic differences, the comparison of 
the results between sBzKs and UV-selected galax- 
ies indicates that the former tend to have slightly 
higher SFRs than UV-selected galaxies with the 
similar stellar masses. This reflects the intrinsic 
difference between the the dust amount in the two 
populations-only sBzKs with higher SFRs may be 
easily detected, due to the larger dust extinction 
than in UV-selected galaxies. 

The positive correlation between SFR and stel- 
lar mass may be attributed to a larger amount 
of gas in more ma ssive galaxies, as suggested by 
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Fig. 9. — The relation between SFRh q and stellar 
mass. Stellar masses are derived from SED fitting 
to BVRi'z'K photometry. The dotted line shows 
the limiting SFR of 87AfQyr~ 1 estimated from the 
lowest detected flux on the assumption that the 
galaxy is at < z >= 1.7 and suffer a typical color 
excess of 0.38. 



5.24. Specific SFR 

We derived specific SFRs (SSFRs) from SFR Hct s 
and stellar Masses. SSFR is an indicator of the 
time scale of star formation activity. Figure [TU1 
shows SSFR as a function of stellar mass. A neg- 
ative relation is seen between stellar mass and 
SSFR, which is co nsistent with previous results 
for z ~ 2 galaxies (|Reddv et al.ll2006at I Erb et al 



2006bl ). This negative relation implies that a 
larger growth of stellar mass o ccurs in les s massiv e 
galaxies. In comparison with lErb et ah (l2006bT) . 
it should be noted that we apply a different IMF 
and dust correction as described in th e previous 
section. In Figure [TU1 the results of I Erb et al 



Daddi et all (|2007l ) 



(|2006bl ) are converted into those with the same 
manner as ours. 

We must consider the effect of flux-limited bias 
on this correlation, as discussed in the previous 
section. However, it seems that there is a genuine 
negative correlation even if the range of validity is 
limited to stellar masses larger than 1 x 1O 1O M0. 
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Fig. 10. — The relation between the specific 
SFRhq and the stellar mass. Filled circles are our 



results, open triangles show those of iReddv et al 



(|2006a[ ) for sBzKs in their UV-selected sample, 
and crosses are those of lErb et all (|2006bl) for UV- 
selected sample. All stellar masses and SSFRs are 
converted into those derived with Salpeter IMF 
and the same dust correction as we apply. The 
dotted line shows the track of galaxies with a con- 
stant SFR of STMoyr" 1 . 



The steep relation between stellar mass and SSFR 
results from the shallow correlation between stel- 
lar mass and SFR, which shows that the increase 
in SFR is not as rapid as that in stellar mass. If 
we miss galaxies with low SFRs of < 87M yr _1 
as discussed in section [5.2.31 that is, the correla- 
tion between stellar mass and SFR is steeper, the 
relation with SSFR would be shallower than that 
in Figure [TOl 

5.3. Equivalent Width 

The rest-frame equivalent widths of Ha are de- 
rived from Gaussian fitting. Both Ha and contin- 
uum fluxes have bee n corrected for du st extinc- 
tion. As discussed in Erb et al. I (|2006bl) . one can 
infer the star formation history of galaxies from 
the equivalent width. While the flux of Ha reflects 
the current star forming activity, the continuum 



Fig. 11. — The relation between specific star 
formation rate and Ha equivalent width in rest 
frame. Rest-frame equivalent widths of Ha are 
derived from line and continuum fluxes after cor- 
rection for dust extinction. 



flux emitted from stars reflects the past star for- 
mation. That is to say, the equivalent width rep- 
resents how actively a galaxy is forming stars com- 
pared with the past. Thus, the equivalent width 
is a similar indicator to SSFR. Figure QT] shows 
the equivalent width as a function of the SSFR, 
indicating a mild correlation. This supports the 
result of the previous section. Less massive galax- 
ies are more actively forming stars, which results 
in a larger growth of stellar mass. 

5.4. Metallicity 

Average metallicity is one of the most impor- 
tant properties of a galaxy. Since gas and stellar 
metallicities are related to the past star forma- 
tion, they are crucial for understanding the star 
formation history and galaxy evolution. It is well- 
known that metallicity is correla ted with stellar 
mass i n the local univers e (e.g 



Trcmonti et al 



20041 ). lErb et all (|2006af ) have reported a simi- 



lar correlation to the local one for UV-selected 
galaxies at z ~ 2. However, their mass-metallicity 
(M-Z) relation at z ~ 2 is shifted to lower metal- 
licity by ~ 0.3 dex. This M-Z relation at z ~ 2 
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obtained by lErb et all (|2006af) is widely used as 
representative of z ~ 2 star forming galaxies 
in comparison w ith both observational and the- 
oretical studies (iLiu et all 20081: Maiolino et al 



2008 : iBrooks et al.l bOoT iDe Dossi et all 120071: 



Finlator fc Da ve 2008). However, because the UV 
selection misses star forming galaxies with a large 
dust extinction, we cannot be con fident that the 
M-Z relation of lErb et aD (|2006al ) applies to all 
star-forming galaxies at z ~ 2. Here, we deter- 
mine a M-Z relation for sBzKs. Our sBzK sam- 
ple, which approximately corresponds to a stellar 
mass-limited sample, should be more suitable for 
deriving the M-Z relation at z ~ 2. 

It is well known tha t electron temperature 
reflects gas metallicit y feewlev fc Dopital 2002 ; 



Kobulnickv fc Kewlevl 12004 lErb et alj l2006al ). 



Because auroral lines, which are used to derive 
electron temperature, are weak, and observed in 
only galaxies with low metallicities, it is very dif- 
ficult even in local galaxies to estimate the metal- 
licity of galaxies from them. We can only use the 
ratio of strong emission lines which are emitted 
from different ionization levels to derive metallic- 
ity for high-z galaxies. At present, there are var- 
ious metallicity diagnostics, such as the R23 and 
N2 methods, to estimate gas-phase oxygen abun- 
dance, 12+log(0/H). One problem with these di- 
agnostics is a disagreement between the resulting 
oxygen abundances inferred from the different di- 
agnostics. Even if the same diagnostic is used, 
some different calibrations give a large difference 
in the resulting metallicities (Kewley fc Dopita 



20021 iKobulnickv fc Kewlevl 12004 ). The same di- 
agnostics and calibration should be used to prop- 
erly compare with other results. 

We derive 12+log(Q /H) using the N2 in- 
dex m ethod proposed bv lStorchi-Bergmann et al 



(jl994h . where N2 is defined as the flux ratio of 
[N II](A6584)/Ha; 



N2 = log 



[NII](A6584) 
Ha ' 



(3) 



The gas-phase oxygen abundan ces are derived 
with t he following relation given in lPettini fc Page! 
(12004b : 



12 + log(0/H) = 8.90 + 0.57 x N2. (4) 

The reasons we use this relation are that it requires 
only two lines of Ha and [N II], and the result 



can be compared directly with that of lErb et al 
(|2006aj ). where the same relation was used to de- 
rive the gas abundance. Also, it should be noted 
that the lines of Ha and [N II] are close enough 
that the difference in dust extinction is negligible. 
The la dispersion i n the abundance deriv ed from 
the relation is 0.18 (|Pettini fc Pagelll2004l ). 

Figure [12] shows the resulting relation be- 
tween stellar mass and metallicity. The broken 
line shows the abundance for galaxies with the 
flux ratio of [N II]/Ha= 0.63, where N2 saturates 
(jKewlev fc Dopital 120021 ) . suggesting that the lin- 
ear relation between the oxygen abundance and 
N2 index is acceptable up to the abundance of 
~ 8.79. As discussed in section^ the three sBzKs 
with a large abundance, ID21193, ID22424, and 
ID31340, may harbor AGNs, because they have 
flux ratios of [N II]/Ha larger than 0.63. How- 
ever, their SEDs suggest that their continuum is 
dominated by stellar emission, and that the AGN 
contribution is small. Even if the three AGN can- 
didates are excluded from our sample, our conclu- 
sions in the following discussion are unchanged. 

We derive upper limits of metallicity for sBzKs 
without a secure detection of [N II] , assuming that 
the peak flux is 3 times sky noise, and that the line 
width is the same as that of Ha. These galaxies are 
shown by crosses in the figure. Also, we estimate 
mean metallicity of the five sBzKs without [N II] 
detection. We derive it from the stacked spectrum 
using the equations (3) and (4), finding that the 
mean abundance is 8.55. This is shown by the 
double circle in the figure, whose stellar mass is 
the average of those of the five sBzKs. 

Figure [T2l shows that more massive sBzKs tend 
to be more metal-ri ch, a trend qualit atively con- 
sistent with that of lErb et all (|2006a[ ). However, 
the metallicities of sBzKs with [N II] detection are 
higher by ~ 0.2 dex than those of UV-selected 
galaxies with the same stellar mass at the similar 
redshift. Considering the uncertainties, a signifi- 
cant shift exists at the 2.7ct level between the two 
M-Z relations. If three AGN candidates are ex- 
cluded from the sample, the confidence level in the 
shift slightly decreases, to the 1.8cr level. 

We again emphasize that the N2 flux ratios 
are converted into t he oxygen abund ances using 



the same relation as lErb et al.l (|2006al ) . We make 



sure that the N2 indices are larger than those of 
Erb et al. ( 2006af ) for galaxies with the given stel- 
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Fig. 12. — Stellar mass vs. metallicity relation. Circles and crosses show our results. Filled circles show 
those of sBzKs with the detection of both Ha and [N II] . Crosses show upper limits of sBzKs without the 
detection of [N II], which are derived using the 3<r upper limit on [N II] flux. The double circle shows 
the mean abundance of the five sBzKs without a detection of [N II] against the mean stellar mass. Open 
squares show UV -selected z ~ 2 galaxies, and open triangles are those of the local ~ 53, 000 SDSS galaxies 



(|Erb et al.ll2006af ). All abundances are estimat ed using the same N 2 index method. The dotted line is the 
solar oxygen abundance of 12+log(0/H)=8.66 (jAsplund et al.ll2004 ). The broken line sh ows the abundance 
for galaxies with the flux ratio of [N II]/Ha=0.63. Stellar masses of the I Erb et alj (|2006a[) are converted into 
those derived with Salpeter IMF. 



lar masses, suggesting that the cause of the differ- 
ence in M-Z relation is not a systematic difference 
in our abundance estimate methodology. 



Pettini k. Pagell (|2004l ) have found a relation 



between gas oxygen abundance and 03N2 index, 
which is defined as 03N2 = log[([0 III] (A5007)/H/3) 
/ ([N II](A6584)/Ha)]. We check the abundance 
with the 03N2 diagnostics for three sBzKs with 
the necessary lines observed. The abundances 
from the 03N2 ratios are lower than those from 
the N2 ratios by ~ 0.13 dex. This is consistent 



with Erb et al.l (|2006af) . who found that the abun- 
dances from 03N2 are 0.17 dex lower on average 
than those from N2. In addition, it is not likely 
that the systematic difference in estimates of stel- 
lar masses accounts for the diff erence in the M-Z 
relation. The stellar masses of lErb et alJ (|2006af ) 



are converted into those derived with Salpeter 
IMF. As discussed in section ISTTl the systematic 
difference of stellar masses should be less than a 



factor of 2. This indicates that there is no way 
that it is the cause of the difference in M-Z rela- 
tion. 

We consider the possibility that our M-Z rela- 
tion suffers from a selection bias that only sBzKs 
with high metallicitics arc plotted, since metallic- 
ities are properly measured only for sBzKs with 
both Ha and [N II] lines detected. The stacking 
analysis for sBzKs without [N II] detection sug- 
gests that there are also sBzKs with lower metal- 
licities (Figure [12]), supporting the idea that se- 
lection bias cannot be completely ruled out. This 
may imply that the mean metallicities of sBzKs 
with given stellar masses are slightly lower than 
those plotted in Figure [12] However, the mean 
metallicity of only sBzKs without [N II] detection 
is similar to that of UV-selected galaxies, imply- 
ing that sBzKs would still be more metal rich on 
average than UV-selected galaxies. This fact sug- 
gests that it is likely that the UV selection misses 
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star forming galaxies at z ~ 2 with comparatively 
higher metallicities, and that metallicities of star 
forming galaxies at z ~ 2 have a large dispersion. 

Therefore, we conclude that the oxygen abun- 
dance of sBzKs is greater than that of UV-selected 
z r~j 2 galaxies. As described in section 15 . 21 
-RT-selected star forming galaxies at z ~ 2 have 
slightly higher SFRs but similar to UV-selected 
galaxies, if the same correction for dust extinction 
is applied, while their metallicities are significantly 
different. This fact may imply that the two star 
forming galaxy populations at z ~ 2 pass through 
different star forming histories before this epoch. 

This difference in the M-Z relation between 
sBzKs and UV-selected galaxies may be related 
to the difference in color of the two populations. 
Figure [T3l shows R — K color as a function of K for 
our sB zK sample and the UV-selected galaxy sam- 
ple of lErb et"aH ( 2006al ). This figure shows that 
the sBzK spectroscopic sample galaxies are redder 
than the UV-selected galaxies. Indeed, the sBzK 
ID31569, which is as blue as UV-selec ted galax- 
i es, is located on the M-Z relation of Erb et all 
The difference in R — K color may imply 
that galaxies with different colors have different 
ages and/or dust contents, since the 4000A break 
falls between the R and -ftf-bands. This fact in- 
dicates that star forming galaxies at z ~ 2 with 
redder colors have larger metal abundances , and 
that the M-Z relation of lErb et all (|2006al ) may 
not represent the entire population of z ~ 2 galax- 
ies. 

Another possibility is that the difference in red- 
shifts of galaxies between the two samples ac- 
counts for the difference in the M-Z relation. 
While the median redshift of galaxies in our sam- 
ple is 1.7, the me an redshift is (z) = 2.2 in the 
Erb et al.l (|2006af ) sample. Because the galaxies 
in our sample are on average at a lower redshift, 
the fact that our sBzKs have higher met a llicities 
than UV-selected galaxies of lErb et al.l (|2006af ) 
may plausibly be attributed to cosmic evolution. 
However, this interpretation would require a sur- 
prisingly rapid evolution of the M-Z relation over 
800 Myr. 

We plot results of star forming galaxies at z ~ 
from SPSS data in F igure [T2l (ITremonti et al 



necessarily descendants of sBzKs. Our M-Z rela- 
tion seems to be simil ar to that of SDSS galaxies 
(jTremonti et al.1 120041 ). but we may not necessar- 
ily claim that the metallicities of sBzKs are the 
same as those of the local SDSS galaxies. There 
is the possibility that the true metallicities of the 
SDSS galaxies are higher than those derived by 
the N2 method, if the saturation of the N2 index 
at high metallicity is taken into account. Thus, 
the oxygen abundances of sBzKs can be as high 
as those of the SDSS galaxies. Then, our result 
that the metallicities of sBzKs are close to solar 
abundance is consistent with the assumption of 
the SEDs with solar abundance in SED fitting. 

It is also found t hat there are M-Z co rrelations 
at other redshifts. Savaglio et al. (|2005l) have re- 
ported the relation for NIR-selected galaxies at 
z ~ 0.7, finding that the metallicities of galax- 
ies with a given stellar mass have a dispersion 
of ~ 0.2 dex around the mean abundances. The 
correlation with large scatter is s imilar to what 
we fo und for galaxies at z ~ 2. iMaiolino et al. 
(|20081 ) have found the relation shifted to much 



2004 lErb et alj l2006al) . It should be noted that 
we cannot directly discuss redshift evolution of the 
M-Z relation, since these SDSS galaxies are not 



lower metallicity for 9 Lyman Break Galaxies at 
z ~ 3.5, suggesting that the metallicity of galax- 
ies decreases with increasing redshift. However, 
since metallicity measurements for high- z galax- 
ies are very limited, especially at z > 2, the M-Z 
relations of high- z galaxies obtained so far may 
have significant statistical and systematic uncer- 
tainties. Therefore, in future, it is essential to de- 
rive more reliable M-Z relations at high redshifts 
from a large, unbiased sample. 

6. Summary 

We have carried out NIR spectroscopy of 44 
BzK galaxies with Kab ^ 23 in the Subaru Deep 
Field with MOIRCS on the Subaru Telescope. The 
z J and HK grism spectra cover wavelengths from 
0.9 to 2.3 //m at a resolution of i?=500. The 
44 spectroscopic targets are selected from a large 
sample of BzKs in 180 arcmin 2 of the SDF. Among 
the 44 targets, 40 are classified as sBzKs and 4 
as pBzKs. Half of the 44 targets are fainter than 
Kab = 22. We detect Ha and some other emission 
lines in 15 sBzKs, nine of which have Kab > 22. 

Spectroscopic redshifts are derived from Ha, 
and fluxes of emission lines arc measured. Stellar 
mass and dust extinction are also obtained from 
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• our sample 

a Erb + 06 sample 

• SDF sBzK sample- 




Fig. 13. — The color magnitude diagram of R — K 
vs. K for our sBzK spectroscopic sample (filled 
circles) , the whole sBzK sample in the SDF (dots) 
and the NIRSPEC sample of lErb et al.l (|2006cj i 
(open triangles). 



SED-fitting to multiwavelength imaging data of 
BVRi' ' z'K. We have then obtained the following 
results. 

The spectroscopic redshift distribution of the 
sBzKs does not peak at z ~ 2, and objects at 
z < 1.4 are also selected. These lower- z objects are 
near the selection boundary in the BzK diagram. 
On the other hand, no z > 2 objects are detected. 
This result suggests that a fraction of star forming 
galaxies at low redshifts can meet the sBzK crite- 
rion. However, the lack of objects at z > 2 may 
be due to a flux-limited bias. 

The Ha luminosities are converted into SFRs. 
It is found that sBzKs with higher stellar masses 
have larger SFRs. The correlation between SFR 
and stellar mass may be attributed to a larger 
amount of gas for more massive galaxies. Divid- 
ing SFR by stellar mass, specific SFRs are also de- 
rived. A negative relation is seen between stellar 
mass and SSFR, which is consistent with previ- 
ous results for z ~ 2 galaxies. This negative rela- 
tion implies that a proportionately larger growth 
of stellar mass occurs in less massive galaxies. 



Gas-phase oxygen abundances, 12+log(0/H), 
are derived from the ratio of [N II](A6584) to Ha 
using the N2 index method. We have found a cor- 
relation between stellar mass and oxygen abun- 
dances, in the sense that more massive sBzKs tend 
to be more metal-rich, which is qualitatively con- 
sistent with the relation for UV-selected z ~ 2 
galaxies. However, the metallicities of sBzKs are 
on average higher by ~ 0.2 dex than that of UV- 
selected galaxies at similar redshifts, and are close 
to those of local galaxies. Compared with the un- 
certainties, this difference in metallicity between 
sBzKs and UV-selected galaxies is significant, and 
may be due to the color difference between the two 
populations. Our sBzKs have redder R — K colors 
than the UV-selected galaxies. Further, an sBzK 
with a blue color similar to the UV-selected galax- 
ies has a lower abundance, comparable to those of 
the UV-selected galaxies. These facts suggest that 
there is a color-dependence in the oxygen abun- 
dance in galaxies at z ~ 2, which may imply that 
galaxies with different ages and/or dust amounts 
have different metallicities. 
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Table 1 

Photometric properties of 15 sBzKs with emission lines detected. 





Magnitude 


Color 


SED Fitting 


ID 


K 


B-z 


z-K 


stellar mass a 


E(B - V) 


21193 


21.3 


1.69 


1.89 


3.8 


0.51 


21360 


21.6 


1.61 


1.53 


1.5 


0.43 


21667 


22.1 


1.24 


2.15 


5.5 


0.38 


21973 


22.4 


1.09 


1.53 


1.9 


0.31 


22424 


22.2 


1.72 


2.65 


15 


0.45 


22602 


22.2 


1.23 


1.98 


5.1 


0.35 


25528 


22.3 


1.21 


1.49 


1.2 


0.41 


26411 


21.7 


1.47 


2.38 


5.4 


0.44 


26821 


21.8 


1.10 


1.24 


2.1 


0.38 


26840 


23.0 


0.73 


0.89 


0.61 


0.29 


27677 


22.4 


1.37 


2.57 


3.6 


0.45 


27815 


22.9 


1.06 


1.45 


1.4 


0.36 


28761 


21.9 


1.04 


1.02 


1.5 


0.35 


31340 


21.6 


1.46 


2.53 


11 


0.41 


31569 


22.0 


0.78 


0.74 


0.61 


0.23 



a The unit is 1O 1O M . 
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Table 2 

Spectroscopic properties of 15 sBzKs with emission lines detected. a 



ID 


Detected Lines b 


redshift 


Fh„ c 


T d 


SFRhq, 


EW Ha f 


12+log(0/H) g 


21193 


Hq,[N II] 


1.479 


10.9±1.77 


1.5±0.24 


4.2±0.68 


1.8±0.56 


8.781°;?* 


21360 


Ha 


1.174 


4.16±2.25 


0.32±0.18 


0.49±0.27 


0.86±0.53 


<8.84 


21667 


Ha 


1.750 


11.1±2.26 


2.3±0.47 


2.5±0.51 


4.5±1.8 


<8.68 


21973 


Ha,[0 III] 


1.504 


18.3±3.58 


2.6±0.51 


1.8±0.36 


8.5±8.4 


<8.75 


22424 


Ha,[N II] 


2.018 


7.25±1.24 


2.1±0.37 


3.9±0.66 


7.0±6.9 


s sn+ 007 


22602 


Ha 


1.872 


7.72±2.35 


1.9±0.58 


1.7±0.52 


1.1±0.46 


<8.88 


25528 


Ha,[N II] 


1.506 


9.99±2.00 


1.4±0.29 


2.0±0.40 


3.3±1.4 


8.94t°;°* 


26411 


Ha,[N II] 


1.498 


14.0±2.67 


2.0±0.38 


3.3±0.63 


2.7±1.2 


8.70tg : » 


26821 


Ha,H/3,[0 III],[N II] 


2.044 


22.2±1.43 


6.8±0.44 


7.3±0.47 


5.9±1.2 


8.66±°:°? 


26840 


Ha,H/?,[0 II], [O III],[N II] 


2.012 


8.98±1.37 


2.6±0.40 


1.5±0.23 


6.6±4.6 


o ro+0.09 
O-O^_0.14 


27677 


Ha 


1.710 


4.53±1.39 


0.89±0.27 


1.6±0.50 


2.8±1.6 


<8.88 


27815 


Ha,[N II] 


1.796 


7.75±1.01 


1.7±0.22 


1.7±0.22 


5.5±4.0 


8.601°;^ 


28761 


Ha,[N II] 


1.818 


10.3±4.90 


2.3±1.1 


2.1±1.0 


1.8±0.91 


8.68t°:Jl 


31340 


Ha,[N II] 


1.719 


7.10±1.39 


1.4±0.28 


2.0±0.39 


2.5±0.73 


o tji+0.08 


31569 


Ha,H/3,[0 III],[N II] 


1.250 


33.0±2.05 


3.0±0.19 


1.1±0.07 


9.9±5.5 


o OO+0.12 
°- /8 -0.25 



a For the two sBzKs classified as an AGN and excluded from the discussions in §5, the SFR, equivalent width, and 
oxygen abundance derived by the same manners as the other sBzKs are also listed. 

b [N II] means [N II](A6584). 

c Observed flux. The unit is 10~ 17 ergs s _1 cm~ 2 . 
d Observed luminosity. The unit is 10 42 ergs s~\ 
°The unit is 10 2 M Q yr _1 . Dust-corrected. 
f The unit is 10 2 A. Rest-frame and dust-corrected. 
g 12+log(0/H) is derived from N2 index method. 
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